Quantum limitations on superluminal propagation 
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Unstable systems such as media with inverted atomic pop- 
ulation have been shown to allow the propagation of ana- 
lytic wavepackets with group velocity faster than that of light, 
without violating causality. We illuminate the important role 
played by unstable modes in this propagation, and show that 
the quantum fluctuations of these modes, and their unitary 
time evolution, impose severe restrictions on the observation 
of superluminal phenomena. 

A famous consequence of Einstein's special theory of 
relativity is the principle of causality: signals cannot 
travel faster than light. Nevertheless, it is known for 
some time, that under certain conditions the group ve- 
locity of an electromagnetic wavepacket can be arbitrar- 
ily large and its energy positive, but yet no conflict arises 
with causality: the signal velocity remains always smaller 
than the velocity of light in vacuum. At least two classes 
of such models allowing "causal" superluminal behavior, 
are known: the first is closely related to tunneling phe- 
nomena, the second arises for unstable systems of atoms 
under an inverted population condition [jj]. 

While in the first case the "superluminal" transmission 
of waves or particles through a barrier is exponentially 
suppressed, the second case shows much more dramatic 
superluminal-like effects. Wavepackets travel with a su- 
perluminal group velocity for unlimited distance, with 
negligible attenuation and dispersion. The second phe- 
nomenon arises due to instabilities in the initial state 
of the radiating system, for example, a scalar field ini- 
tially in a "false vacuum" state Q|. Recently, an op- 
tical experiment studying superluminal group velocities 
for waves in a medium with inverted population was sug- 
gested [^|-||. While the discussion was limited so far to 
the semi-classical limit, it was further suggested ||, that 
with a full fledged quantum-mechanical treatment one 
would possibly obtain stable tachyonic-like quasiparticle 
excitations in the inverted medium. 

In this Letter we examine classical and quantum as- 
pects of superluminal group velocities in unstable sys- 
tems. Classically, a peak of an analytic wavepacket may 
travel from point A to point B faster than the speed of 
light, since its shape at B can be fully reconstructed from 
the part of the wave packet which is causally connected 
to the point B. We show that quantum mechanically, 
analyticity of the wave packet is not the only condition 



needed for such a reconstruction. Rather, the part of the 
wave packet which is causally connected to the point B 
must contain many photons or large enough energy in 
the unstable modes. This condition strongly suppresses 
superluminal effects in the limit that the system con- 
tains only a few photons. The exponential suppression 
of superluminal effects characterizes both the tunneling 
phenomena and the present case of unstable systems. 

To begin with we introduce a simple model which clas- 
sically exhibits tachyonlike motion. Consider a real scalar 
field <p(x,t) in one spatial dimension, under the Hamil- 
tonian 
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Here, n(x, t) = <p{x,t) is the field conjugate to ip(x,t), 
A > and the speed of light is put equal to one. This 
Hamiltonian describes the continuum limit of coupled 
pendula. In the stable ground state <p(x) = tt. We will 
be interested, however, in the dynamics of ip near the 
metastable state <p(x) = 0. For any finite time interval of 
interest, we may take A sufficiently small such that the 
potential term may be expanded, cosa/A^ « 1-|^ 2 . 
From now on we restrict ourselves to this expansion. 
Then, the equation of motion for tp becomes, 
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In term of the cigenmodes the solution is 
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where reality conditions imply (p~k — ipt, n-k = tu.- 
The time evolution of the modes is given by 



<Pk(i) = <fiko cosw K i + — sinw K i 
7Tfe(t) = 7r fc coscjfci - ojkfok sinwfci 



where 
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Notice that for |fc| > m, u>k is real, and ipk,^k oscil- 
late in time. In the range \k\ < m LUk is imaginary, and 
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ifk(t) and TTk(t) are exponentially diverging. The latter 
modes are analogous to spontaneous emission in the op- 
tical model of an inverted medium of two-level systems. 
We will henceforth refer to the free oscillatory modes and 
diverging modes as the "normal" and "unstable" modes, 
respectively. It is only near the point ip ~ tt ~ that 
these modes exist as linearly independent solutions. As 
tp grows sufficiently the instability is damped by the non- 
linear A-term in ([!]), but for the time of interest to us, 
the latter can be neglected. 

We now turn to examine several features of the classical 
equation of motion (||). First, we examine the propaga- 
tion of a wavepacket given by 

<p(x,t) = f dk(g (k)e ikx - iUkt + h.c.) (7) 

J\k\>m V ' 

and n(x, t) = dt<p(x,t). We take go(k) to be non zero 
only in the range of normal modes m < \k\ < k max , and 
centered around k = ko, with a width Ak <C ko — m. The 
spatial width of this wavepacket is 1/Afc, and we assume 
that at t — it is spatially centered very far to the left 
of the origin, at Jo < 0. 

This wavepacket propagates with a group velocity, v g , 
and a phase velocity, v p , given by 

v 9 = — = m 9 ( 8 ) 
v P ^/kl - to 2 

Since v g > 1, the motion of the center of the wave packet 
(group velocity) is superluminal (tachyonlike) , while the 
phase velocity is always subluminal. As long as the dis- 
persion is negligible, <p(x, t) is (up to a time dependent 
phase factor) just <p(x—v g t, 0), i.e., the initial wavepacket 
moving at velocity v g . 

Second, we note that despite the superluminal group 
velocity, causality is maintained ^,||. This is best seen 
through the Green's functions of Eq. @. In terms of 
the homogeneous Green functions: (p(x,t) = J dx'[G(x — 
x' ,t)ip (x')+G(x—x' ,t)n (x')]. However, G and G vanish 
outside the light cone, i.e., for x — x' > t, therefore the 
value of <p at point x at time t is fully determined by its 
value at points x' that are causally connected to x. 

For further insight into this point, it is instructive to 
examine what happens if the initial wavepacket is trun- 
cated at x — 0, i.e., ipo(x) and ttq(x) are replaced by 
(Pq = 0(x)ipo(x) and ttq = 9(x)tto(x), respectively. Here 
9 is a smoothed step function. The length scale over 
which it is smoothed is assumed smaller than all other 
length scales in the problem, and is kept finite just to 
avoid an infinite d x <p. Note that since ipo(x) is centered 
at Xq <C 0, the truncated wavepacket, <Pq(x) consti- 
tutes only the tail of the original one, and the energy 
stored in the truncated wavepacket is just a small frac- 
tion of that stored in the original one. The momen- 
tum representation of the truncated wavepacket tpo(x) is 



9o(k) = / dk' k _fci +iri gc(k'), with 77 an infinitesimal num- 
ber. The smoothing of the step function is accounted for 
by introducing an upper cut-off to the k integral. While 
the original wavepacket was solely composed of normal 
modes, the truncated one includes also unstable compo- 
nents. 

The time evolution of the truncated wavepacket can 
be calculated using the exact expression for the Green 
function ^,||||. There are three regimes: (a.) For points 
x < —t, the causality of the Green function dictates that 
the value of ipx{x,t) is not affected by the existence of 
the truncated wavepacket at x > 0. Thus, the field is 
zero in this range, (b.) For — t < x < t, however, 
the time evolution of the truncated wavepacket is very 
different from that of the original one: it is exponen- 
tially growing due to the contributions of unstable com- 
ponents. Since the amplitude of the normal modes oscil- 
lates in time while the amplitude of the unstable modes 
grows exponentially, for times t 3> 1/m the wavepacket 
tp T (x,t) is predominantly composed of unstable modes 
of very small wavevectors, with an exponentially small 
contribution of normal modes, (c.) Finally, for points 
x > t, the causality of the Green function dictates that 
the value of tpx(x,t) does not depend on whether or not 
it is truncated at the origin at t = 0. It is given, up 
to a time dependent phase factor, by ip(x — v g t, 0). The 
propagation of the superluminal peak is not affected by 
the truncation, and, for long enough times, one finds at 
Xo+v g t a wavepacket of width 1 / Ak, in which the field ip 
oscillates with a characteristic wavevector of kg. Thus, at 
long times the field <p T {x, t) shows an interesting behav- 
ior: it is predominantly composed of very small wavevec- 
tors \k\ < m with just an exponentially small contribu- 
tion of wavevectors |fc| > to. However, in a region of 
spatial width 1/ Ak around Xq + v g t it oscillates at the 
large wavevector ko > to. The region 1/ Ak can be made 
arbitrarily large, if t is taken long enough. These short- 
wavelength oscillations of a long-wavelength wavepacket 
are very similar to the superfast Fourier oscillations dis- 
cussed in [0 [o) . 

The above discussion demonstrates that in order to 
maintain consistency with causality, the mechanism 
which gives rise to superluminal group velocities has to 
rely on the local information stored in the tail. As long 
as the true information that the wave has been truncated 
has not arrived, the local amplification extrapolates the 
peak as if this truncation does not exist. Classically, 
the possibility of extrapolating the full wave from an in- 
finitesimally small tail is, although surprising, possible. 
As discussed in ||, for an analytic function, the wave 
can be reconstructed by means of a local Taylor expan- 
sion, and the tachyonlike propagation can be viewed as 
an analytic continuation of the tail. 

As we now show, the extrapolation of the full propa- 
gated wavepacket from the truncated tail is made possi- 
ble by the unstable modes, and this is true independent 
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of the details of the model. Since g${k) is the Fourier 
transform of a small tail, it is very small. The propa- 
gated wavepacket, after truncation, is 



if T (x, t) = / dk ( e 



ikx — iuikt T 
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However while on the right hand side, the amplitudes 
§q (k), originating from the truncated tail, are very small, 
on the left side the magnitude of the tachyonic peak is 
not small, since the peak is fully reconstructed from the 
tail. Since the amplitude of the stable modes, ipr(k), for 
|fc| > m, is constant in time, this amplification can arise 
only from the contribution of unstable modes growing 
exponentially with time. It is then essential to have un- 
stable modes for a full reconstruction of a superluminal 
peak from a truncated tail. 

We now proceed to the quantum mechanical analog of 
the classical superluminal system. Can a similar mech- 
anism of local amplification give rise to superluminal 
group velocities? If so, it is conceivable that by means of 
a process, analogous to the classical analytical continu- 
ation, this local information may be used to reconstruct 
the wavepacket's peak. 

The causality of the Hamiltonian ([!]) is manifested in 
the quantum case by the statement that causally discon- 
nected local observables commute. Again, this causality 
motivates us to study the quantum analog of the classi- 
cal truncated wavepacket. In this study, detailed below, 
we define an analog of a "false vacuum" state for this 
model, on top of which we define a quantum state de- 
scribing a truncated wavepacket. Then, we show that if 
the truncated wavepacket is too small (in a sense to be 
explained below), the initial quantum state describing it 
is not orthogonal to the false vacuum state. In that case, 
the vacuum quantum fluctuations in the fields tp, ~n are, at 
t — 0, larger than the amplitude of the wavepacket, and 
one cannot distinguish between the false vacuum state 
and the wavepacket state. As we saw above, classically 
the small truncated wavepacket becomes exponentially 
large in time, and one may naively hope that quantum 
mechanically it eventually becomes larger than the vac- 
uum fluctuations. However, this cannot be the case, due 
to the unitary quantum mechanical time evolution. By 
the principle of unitarity, the scalar product of two states, 
here the false vacuum and truncated wavepacket states, 
is constant in time. Thus, if at t = the truncated 
wavepacket is not orthogonal to the false vacuum, it can- 
not be distinguished from it at all later times. The ex- 
ponential growth of the tail is masked by an exponential 
growth in the quantum fluctuations. 

Turning to a more concrete discussion, we promote the 
functions <p(x, t), tt(x, t) to quantum operators. From the 
standard canonical commutation relations for the fields, 
[tp(x, t), ir(x' , t)] — iS(x — x'), it follows as usual that 
for the normal modes: [^k,7r|,] = ""fc'] — i8(k — k'). 



However, for the unstable modes the the system does not 
admit the usual Fock space structure with an ordinary 
particle interpretation. For these modes (|fc| < m), the 

Hamiltonian, J^k 7r I 7r fc — ( m2 ~ k 2 )ip\.ipk, is unbounded 
from bellow, and has a continuous and unbounded spec- 
trum of energies for each wave number k. Consequently, 
these modes cannot be put in a ground state, and the 
fluctuations in the fields (pk , i^k inevitably grow exponen- 
tially with time. 

For the normal modes \k\ > m the vacuum state is 
obviously the ground state of an harmonic oscillator of 
frequency \Jk 2 — to 2 . The situation is less obvious for 
the unstable modes. On one hand, we would like the 
quantum fluctuations in the fields (pk , i^k to be as small 
as possible, but, on the other hand, stationary states of 
the unstable modes necessarily have infinite fluctuations. 
Any state of finite fluctuations is non-stationary, with 
the fluctuations growing exponentially in time. Thus, we 
choose the initial state of the unstable modes to be a 
non-stationary one, in which the exponential growth of 
the field fluctuations is slowest. This state is the ground 
state of an harmonic oscillator of frequency \Jm 2 — k 2 . 
The vacuum state of the entire system is then a direct 
product of the vacuum state for each mode k. 

Next, we should find the quantum analog of the wave 
packet we analyzed in the classical case. Again, causal- 
ity dictates that a measurement of the field ip at a given 
point cannot be affected by its values in causally discon- 
nected points. Thus, to analyze the quantum analog of 
the classical super-luminal propagation, we need to con- 
struct the quantum analog of the truncated wavepacket 
lPq(x). We do so by shifting each mode ifk from its vac- 
uum state, in which (ipk) = 0, to a coherent state, in 
which ((pk) = go(k): 



l*o> = exp 
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where |vac) is the vacuum state. Classically, as we saw 
above, the truncated wavepacket included all the infor- 
mation needed to reconstruct the superluminal propaga- 
tion of tpo(x). Quantum mechanically, however, this is 
true only if the state (jJ^j is orthogonal to the vacuum 
state. If that is not the case, one cannot distinguish be- 
tween the time evolution of the vacuum state and that of 

©■ 

The scalar product (vac|\I/o) can be easily calculated. 
It is, 



(vac I* 
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It is useful to distinguish between the normal (|fc| > m) 
and unstable (|fc| < to) contributions to the /c-integral 
in (|ll|). The physical interpretation of the normal con- 
tributions is rather clear: g$ (k) is the amplitude of the 
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oscillations of the mode k, and oj(k)\(gQ (k)) 2 is the aver- 
age number of photons in that mode. The contribution 
of the unstable modes cannot be discussed in terms of the 
photons. For these modes, the integrand i \u>(k)\ (g^ (k)) 2 
is just the ratio of the kinetic energy stored in the fc-mode 
to Tku(k). The quantum state containing the wavepacket, 
| ^o) > is then orthogonal to the vacuum state only if the 
truncated wavepacket contains many photons (^$> I) of 
the normal modes, or large enough energy in the unsta- 
ble modes. This poses a new condition, which did not 
exist in the classical theory, for the reconstruction of a 
tachyonic peak from a truncated tail: the truncated tail 
cannot be too small. If it is too small, the initial fluctua- 
tions in the field ip, which grow exponentially with time, 
overcome the reconstructed peak. 

To exemplify the last point, let us evaluate the fluc- 
tuations of a local observable. The field at a point has 
singular fluctuations, and one needs to consider smeared 
operators like: P(L,t) = J f(x')ip(L — x',t)dx'. where 
fix'), is non- vanishing only in the spacetime volume AL 
around the point L, and is normalized to / fdx = 1. 
Since (\I/o I I ^o) is identical to the classically smeared 
field, the fluctuations are dominated by (4'o|^ 2 | v I'o) = 
(<P 2 )n + (<p 2 )u where the subscripts N and U, stand for 
the contributions of the normal and unstable modes. The 
contribution of the normal modes is always finite. On 
the other hand, the unstable band yield an exponentially 
growing fluctuation: 



)u^hI [2mT] T>> -X /m h 
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where Iq is zeroth order modified Bessel function. For the 
observation of superluminal propagation two conditions 



have to be satisfied: First v g T ol 



> (where T obs is the 



time at which the wavepacket is observed), i.e., the point 
of observation should be far outside the initial spread of 
the wavepacket. Second, we must distinguish between 
superluminal propagation and propagation at the speed 
of light. This leads to: (v g — l)T b s ^> jh:- These con- 
ditions require that T m i n m 3> 1. Eq. (12) thus implies 
that for the signal amplitude to be larger than the am- 
plitude of the fluctuations at the observation time, the 
signal amplitude should ne exponentially large. 

One may still question the generality of our result. 
For instance we have assumed that initialy the unsta- 
ble modes are in a direct product state. Can a dif- 
ferent initial state be made orthogonal to the vacuum 
more easily? We now argue that the obstacles pointed 
above in the way to realizing quantum mechnical super- 
luminal group velocities cannot be circumvented. Con- 
sider three initial classical tachyonlike wavepackets, de- 
noted by R, L and D, with intial field and momenta 
configurations (tp R , n R ), (ip L = (p R , tt l = —itr), and 
(<Pd = —fR^D = 7!\r)- R and L represent right and 
left moving wave packets with the same initial <p, D is 



right moving with negative (p. Suppose that initially the 
wavepackets are localized around Xq <C 0. Then after 
some time, R and D form peaks (of opposite sign) at 
X 3> while L leaves at X ^> only a small tail. Again, 
we consider the quantum states analogous to the clas- 
sical truncated tails of the three wavepackets, denoted 
by R T ,L T ,D T , where the truncation is done at X = 0. 
Classically, all three wavepackets can be reconstructed 
from their tails. Quantum mechanically, however, it is 
not possible for all three quantum states to be mutually 
orthogonal: R T and L T are orthogonal if the uncertainty 
in the momentum in each state is smaller than the dif- 
ference between the momentum expectation value of the 



two states: 



— (7r£). Similarly, R and 

D are orthogonal if A(pi RD -, <C — However, 
at the tail, the right hand side of these two inequalities 
can be made arbitrarily small, while the left hand side 
cannot be made arbitrarily small, due to the uncertainty 
principle. Thus, if the tails are small enough, the classi- 
cally distinguishable wavepackets R, L and D cannot be 
mapped into three mutually orthogonal quantum states. 

In conclusion, we have shown that classical superlu- 
minal - like effects become incompatible with unitarity 
in the quantum mechanical limit, and are strongly sup- 
pressed. Our argument strongly questions the possiblity 
that these systems may have tachyon-like quasiparticle 
excitations made of a small number of photons |j) . 
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